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dominates in the Regge regime. The graphic representation of the leading trajectory is very similar 
to the diagrams widely used to illustrate meson scattering. 
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In this Letter we study the large N QCD scattering 
amplitudes of colorless composite quark operators given 
as a path integral over Wilson loops [1, 2]. We find the 
most important paths by means of a double saddle point 
method. A geometric representation of these paths is 
very similar to the dual diagrams which for many years 
have been used to illustrate meson scattering amplitudes. 

Previously we have considered [3] the large N (or 
quenched) QCD meson scattering amplitude, using the 
connection between this amplitude and the Wilson loop 
in the case where the loop is well represented by the area 
behavior. This QCD scattering amplitude then turns out 
to be the Regge limit of the usual Veneziano amplitude in 
the cases when the number of external particles is large, 
or when the quark mass is small relative to the QCD-A. 
As discussed in [3] this is a prediction for near forward 
scattering at high energies, that is associated with large 
distances in the relevant functional integral. The goal of 
this Letter is to demonstrate this by explicitly construct- 
ing the saddle point trajectory. 

The relation between the off shell M meson scattering 
amplitude G and the Wilson loop occurs as a functional 
integral [3], 
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in the limit of large M or small quark masses, when the 
spin dependence factorizes. In this expression W{x{(j))) 
is the Wilson loop defined on the boundary curve x{(f>), 
and the function p{(f>) is piecewise constant, defined by 



p{(t>) = for (f>i < (f>< (f>i+i, 



(2) 
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where the momenta qi are related to the external mo- 
menta Pi's by 



Pi = qi-i - qi 



(3) 



thus ensuring conservation of the energy-momentum. In 
particular we note that Eq. (2) implies the very useful 
result 



(4) 



Also, in Eq. (1) the Wilson loop is required to be bounded 
by the closed curve a;^ (</)). 

In the basic relation (1) we now want to insert the 
area behavior W ~ exp(-area), where the area is mea- 
sured in terms of the string tension. If we insert the 
area as the usual square root of the determinant of the 
extrinsic metric and minimize this, the functional path 
integral in (1) becomes impossible to perform because it 
depends on surface variables instead of the paths in (1) 
and also because of the highly non-linear character of the 
integration. We therefore suggested [3] to use Douglas' 
expression for the minimal area [4] , 
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This integral should be minimized with respect to the 
reparametrization function 9{(j)), subject to the require- 
ment 9{(j)i) = (f>i. We should therefore look for a saddle 
point with respect to the reparametrization. 

The expression (5) is known from the derivation of the 
string disk amplitude in the Polyakov formulation, where 
it appeared as the boundary action after integration over 
fiuctuations of an open string with the Dirichlet bound- 
ary conditions at the ends. It is less known that the 
crucial role is then played by the path integral over the 
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boundary metric (the same as the path integral over the 
reparametrizations), which is dominated for large loops 
by the saddle point in hands, reproducing the minimal 
area. However, for small loops W{x{(j))) in Eq. (1) is 
described by QCD perturbation theory, rather than by 
the string disk amplitude, thus leading us to the under- 
standing [5] of QCD string as an effective rather than 
fundamental string. It makes sense only for long strings 
and we shall soon justify by an explicit calculation that 
the double saddle point we consider is indeed associated 
with large loops or long QCD string. 

Douglas' integral (5) is apparently well suited for in- 
sertion in the path integral (1), where it leads to an 
apparent Gaussian integration. However, this point of 
view is an illusion - the non-linearity is still with us, be- 
cause we have to minimize the integral (5) with respect 
to reparametrizations 9{(j)) or the inverse ^{6). As shown 
by Douglas, this leads to the requirement 
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This condition is highly non-linear! 

We now proceed to tackle this problem by a two 

step mechanism: first we find the saddle point in the 
x— integral for any reparametrization (justified because 
we seek a large minimal area), and then we require that 
for this saddle point the reparametrization miist satisfy 
the condition (6). Thus we search for a double saddle 
point. 

From Eq. (1) the saddle point condition for large area 
or, alternatively, large string tension is 

/■27r 

ip^{0) + — d9' xf'ie') log[l - cos{(j){e) - (A(6l'))] = 
Jo 

(7) 

for W ~ exp(— crA), where a is the string tension. This 
equation can easily be solved by noting that the log is a 
Green function with a simple inverse, 
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G-i(^) = -G(0) = 
Inverting (7) we therefore get 
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W = 1^ r de'p'^ie') iog[i - cos(<^(^) - cp{9'))]. 

(9) 

Note that this saddle point x is imaginary, as is to be 
expected from the occurrence of ipx in (1). Next we can 
use Eq. (4) to obtain 

x'^m)) = ^YlPi ios[i - - '^j-)]- (10) 

We note the subtle point that in principle the right-hand 
side depends on the reparametrization, but the left-hand 



side is independent! Thus, the saddle point does not 
care about 0{(j)), i.e. the reparametrization can be any 
function. 

We now have to show that the saddle point (10) actu- 
ally satisfies the non-linear Douglas condition (6). From 
(10) we have 
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Inserting this on the left-hand side of Eq. (6), we obtain 
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The left hand-side can be easily evaluated if we perform 
a partial integration and use Eqs. (8a) and (8b), 

left-hand side of (12) 
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In order to satisfy the Douglas condition (6) we therefore 
need to have 
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for each value of j. Among these equations only M — 3 
are independent because of the projective symmetry, thus 

fixing M — 3 integration variables (j)i 's as a consequence 
of the fact that we arc looking for a double saddle point. 
The remaining three variables can be fixed arbitrarily. 

Instead of the angle (j) ranging between and 27r it is 
often convenient to use 



r = -cot((/)/2), 



(15) 



which ranges from —00 to -t-00. The saddle point solution 
(10) then becomes 



(u(r)) = — V'p^ log \r - ■ 
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Here u{r) is any reparametrization of the boundary 
curve. 

Since the trajectory (16) is the saddle point of the am- 
plitude (1), it can be approximated by the value of the 
integrand at this saddle point. We thus obtain [3] 



^—aA[x]+i J p{r)x{r)dr\ 
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= g(l/27r(T) PiPi log|ri-rj| 



(17) 



which coincides with the saddle point approximation to 
the integral over the Koba-Nielsen variables in the M 
point Veneziano amplitude. 
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This steins from the fact that the Douglas procedure 
leads to a determination of the parameters r^'s in a sim- 
ilar way how we have seen this in Eq. (14) for the corre- 
sponding parameters (/)i's. Proceeding as before we ob- 
tain [6] 
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The meaning of this condition is easy to understand: on 
the right-hand side of Eq. (17) the condition (18) corre- 
sponds to the extremum, generally including the Regge 
limit, and this extremum is thus exactly the minimal sur- 
face entering the Wilson loop in QCD! 

It is instructive to see how the curve (16) looks like 
for the four point function in the center of mass frame. 
Because we deal with the case where both Mandelstam's 
variables s and t are negative, this is associated with a 
scattering process in u-channel. In view of the projec- 
tive symmetry we can set three out of four r^-'s equal 
to certain typical values, keeping in mind that one such 
contour belongs, in fact, to a whole family of contours, 
which are equally important. The projective invariant 
ratio (r2 — ri)(r4 — r^) / {r^ — ri)(r4 — r2) is then fixed by 
Eq. (18) to be t/{t + s). We plot the corresponding con- 
tour in Fig. 1 for s/t = 10 in the units of iypaj^-Ka. The 
contours are quite normal for not too small tj s and look 
pretty much similar to dual diagrams used to illustrate 



the scattering of mesons. For r 



we have 



-•p^A X regularizing factor. (19) 
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Actually, Fig. 1 is obtained without any regularization, 
which is implemented by Mathematica itself. The 
transverse size of the contour shrinks when f/s — s- in 
spite of the fact that the minimal area 



o'5'ini„(C) = a'tlog-, 
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can be large even for small t. For the self-consistency of 
our approach this transverse size should be larger than 
the confinement scale 1/Aqcd ~ 

A very interesting question is as to how we can make 
the saddle point trajectory (16) to be real. This happens 
in the physical domain of s-channel, when s > 0, i < 
and M < 0, so that —i\/u = \fs. We may think about this 
as an analytic continuation in s from negative to positive 
values, keeping |s| —t all the way. The area of such a 
surface will be a complex number because of Minkowski- 
space kinematics but with a negative real part, as is seen 
from Eq. (20), what is essential in our consideration for 
it to dominate in the path integral. 

There is a number of well-known reasons why the am- 
plitude (17) cannot be right in all kinematical regions 
of QCD. For example, for large transverse momenta the 
decay should be power like, in contrast to the dual mod- 
els, and furthermore, QCD has of course no tachyon. In 
this connection it should be remembered that the area 








FIG. 1: Saddle-point trajectory (16) for tja = 0.1: 3D view 
and its projections in 0,1; 0,2 and 1,2 planes, respectively. 
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behaved Wilson loop is quite small, so it may only dom- 
inate in restricted kinematical regions. For large trans- 
verse momenta the amplitude corresponding to Eq. (17) 
decays exponentially, so the power behaved perturbative 
behavior of W is much more important for large trans- 
verse momenta. We therefore expect the area behavior to 
dominate only at large energies in the near forward direc- 
tion. Since the tachyon is a low energy phenomenon [5] 
it is therefore not relevant in our work. 

From Eq. (16) we can obtain the equation for the min- 
imal surface because x has to be harmonic, 

x''{r,y) = —'S2p'i^og\z-rj\, z = r + iy, y>0. 
Tver '—r 

(21) 

This expression is the imaginary part of the analytic func- 
tion 

— ^ppog(z-r,), (22) 

3 

whose real part 

— V arctan (23) 
-KG ^ ^ y ^ ' 

is also harmonic and describes [6] the minimal surface 
bounded by a step function 

i^p^^e(r-r,)e(r,+i-r) (24) 
i 

with being the Heaviside step function. It is of course 
related to the imaginary part by the Cauchy-Riemann 
equations, which in turn are related to a T-duality map 
[7]. 

The induced metrics of the surfaces (21) and (23) dif- 
fer only by the signs. Moreover, if we were go beyond 
the saddle point approximation, their transverse quan- 
tum fluctuations can be similarly taken into account by 
a conformal mapping of the upper half-plane onto a rect- 
angle, associated with the worldsheet parametrization. 
We may expect therefore that we shall obtain for (21), 



like in [6] for (23), the linear Regge trajectory 

"W = ^+"'* (25) 

in the semiclassical approximation for arbitrary dimen- 
sion d > 2. 

The expression (17) leads to straight line Regge tra- 
jectories, when the values of ri's obey Eq. (18). The 
prediction is thus that this should occur at high energies 
and moderate transverse momenta, in accordance with 
experimental results for SU(3). Also our results are valid 
in other dimensions than four, for example in two dimen- 
sions, where 't Hooft [8] found that the mass spectrum for 
large masses is a linear function of an integer quantum 
number. For smaller masses there is a deviation from 
the straight line. In two dimensions the Wilson loop is 
area behaved everywhere, but our calculation only in- 
cludes the saddle point. For smaller areas there would 
be corrections to what we found. 

The fact that the minimal surface (21) is associated 
with the Regge limit of the Koba-Nielsen amplitudes is 
known since early days of string theory [9], though the 
reparametrization of the boundary was not taken into ac- 
count that time. To get it within our approach, it was 
crucial that the path integral (1) is dominated by the 
saddle point trajectory (16) whose shape is not sensitive 
to the choice of the reparametrizing function u(r). It 
is also worth mentioning another minimal surface of the 
form of a helicoid that was associated in Ref. [10] with 
the reggeization of quark-antiquark within the AdS/CFT 
correspondence in a confining background. It will be in- 
teresting to understand how these two are related. How- 
ever, we emphasize once again that our approach is based 
solely on the representation of the Wilson loop via the 
boundary functional (5) and does not use explicitly a 
surface representation. 

As a conclusion we see that the large N meson scat- 
tering amplitude expressed in terms of a path integral 
is dominated by a double saddle point for large energies 
and moderate transverse momenta. The graphic repre- 
sentation of the leading trajectories is similar to the dual 
diagrams used to illustrate meson scattering. 
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